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We determine the statistics of work in isothermal volume changes of a classical ideal gas consisting 
of a single particle. Combining our results with the findings of Lua and Grosberg [J. Chem. Phys. 
B 109, 6805 (2005)] on adiabatic expansions and compressions we then analyze the joint probability 
distribution of heat and work for a microscopic, non-equilibrium Carnot cycle. In the quasi-static 
limit we recover Carnot efficiency, however, combined with non-trivial distributions of work and 
heat. With increasing piston speed the efficiency decreases. The efficiency at maximum power stays 
within recently derived bounds. 
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I. INTRODUCTION 

Thermodynamics has to be modified for systems so 
small that typical changes of their energy are of the order 
of the thermal energy per degree of freedom. To account 
for the correspondingly strong fluctuations of energy, 
work, heat, and entropy, in the emerging field of stochas- 
tic thermodynamics thermodynamic quantities are char- 
acterized by probability distributions [T]-[3] . These distri- 
butions exhibit unexpected symmetries, now subsumed 
under the notion of fluctuation theorems, which have re- 
vealed new insights into the role of large deviations in 
statistical mechanics and the foundations of macroscopic 
thermodynamics. Thanks to recent advances in experi- 
mental techniques these findings are amenable to exper- 
imental verification and contribute to the understanding 
of energy conversion on the nano-scale. A particularly in- 
triguing application concerns the efficiency of biological 
motors [H [5] as well as miniaturized heat engines [B| . 

One of the remarkable features of fluctuation theorems 
is that they hold for (almost) arbitrary deviations from 
thermodynamic equilibrium. In the attempts to intu- 
itively understand this surprising generality case studies 
of exactly solvable model systems have played a decisive 
role [THS]. Given the ubiquitous presence of the ideal 
gas model in statistical physics it is not surprising that 
among these studies those making use of ideal gases P" 
113] are particularly prominent. In 2005 Lua and Gros- 
berg [9] scrutinized the validity of the Jarzynski equal- 
ity [13] for the adiabatic expansion and compression of 
a classical ideal gas consisting of one particle and high- 
lighted the importance of the far tails of the Maxwell 
distribution. Soon afterwards an analogous treatment of 
isothermal volume changes was pursued [12] . however, 
analytical progress was limited to the approximate treat- 
ment of the case of large piston speed and small volume 
change. To characterize the general case only numerical 
simulations were used. 



In the present paper we re-analyze the isothermal ex- 
pansion and compression of a single classical gas particle 
in the framework of stochastic thermodynamics and show 
how some of the obstacles hampering analytical progress 
can be overcome. Analogous to Lua and Grosberg we 
consider the particle being inside a cylinder with a move- 
able piston and make no a-priori assumptions about the 
piston speed such that also strong deviations from equi- 
librium are possible. We determine the distribution of 
work performed by the gas particle as well as the statis- 
tics of the heat exchanged with the reservoir. As a test of 
our findings we check their consistency with the Jarzyn- 
ski equality and the Crooks fluctuation theorem [15] . 

We then consider a cyclic process of Carnot type con- 
sisting of two adiabatic and two isothermal strokes. Com- 
bining the results of Lua and Grosberg for the adiabatic 
parts and our own findings characterizing the isothermal 
changes we determine the joint distribution of work and 
exchanged heats for one cycle of the engine. We show 
that even in the quasi-static limit of slow piston speeds 
fluctuations prevail and that the properties differ from 
those of the classical Carnot process. Finally, with the 
distributions of work and heat at hand, we analyze the 
efficiency of the engine at maximal power and compare 
the result with previous studies [THl [13 and with recently 
derived bounds [THH^S]. 

The paper is organized as follows. In section [llj w e in- 
troduce our model and fix the notation. Section lllll con- 
tains the analysis of an isothermal expansion and com- 
pression and parallels the treatment of Lua and Gros- 



berg for adiabatic processes. In section IV a full cycle 
consisting of two isothermal and two adiabatic strokes is 
considered. We first investigate the quasi-static limit in 
which the piston speed is small and then turn to the gen- 
eral case of arbitrary piston speed including a discussion 
of efficiency at maximum power. Finally, we summarize 
our results in section W\ 



II. THE MODEL 
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The model consists of a one dimensional cylinder of 
length L closed on top by a piston moving at constant 



2 



speed u. The cylinder contains one classical particle with 
mass m = 1 that moves with velocity v. At the piston 
the particle is assumed to be reflected elastically. With 
the piston being much heavier than the particle the ve- 
locity of the particle is changed from v to 2u — v upon a 
reflection whereas the piston speed remains unchanged. 
The work done on the system by one collision of the par- 
ticle at the piston is given by the change of the particles 
kinetic energy 



AW^ = ~2u{v - u). 



(1) 



In case of an adiabatic process [S] the particle is re- 
flected elastically at the bottom of the cylinder as well. 
Consequently there is no heat exchange in this case. To 
describe an isothermal compression or expansion we ther- 
malize the particle in the collision with the bottom of the 
cylinder, i.e. upon reflection we draw a new random ve- 
locity from the distribution [T^ [TB] 



= /3ve ^ l[v > 0] 



(2) 



Here Boltzmann's constant is set to unity, /? denotes the 
inverse temperature and the indicator function ![•] is one 
if the argument is true and zero otherwise. At the bottom 
the particle therefore exchanges heat with the attached 
reservoir but no work is performed. The distribution ([2| 
is chosen such that, for a static piston, u = 0, the equilib- 
rium Maxwell distribution is preserved. For constant /3 
we call the process "isothermal", although, for a moving 
piston, in particular with large w, it is in general impos- 
sible to associate a temperature with the gas particle. 



III. ISOTHERMAL EXPANSION AND 
COMPRESSION 

A. General analysis 

In this section, we discuss an isothermal change from 
the initial cylinder length Lq to the final length Lf = Lq + 
uT during time T with constant speed u of the piston. 
At time t = the system is considered to be in thermal 
equilibrium, i.e., the distribution of initial position xq 
and initial velocity vq of the particle is given by 



Po(a;o,wo 




exp 



-13- 



l[0<xo<Lo]. (3) 



To determine the distribution of the total work W 
transferred during an expansion or compression, we need 
to know the times tj at which the particle hits the bot- 
tom. With the help of ^ we may then determine its 
velocity and, using Q, also the work performed at the 
piston. If the particle leaves the bottom at time tj with 
velocity Vj, it reaches the bottom again after the time 
increment Atj = 2Lj/(vj — 2u). Here Lj = Lq + utj is 
the length of the cylinder at time tj . The time interval 
Atj not only depends on Vj but, through tj, also on all 




Figure 1: Classification of particle trajectories on the basis 
of the initial position xo and the initial speed vq. For either 
Vq < «Q^' or iiQ > iIq^' the particle reaches the bottom, x — 0, 
of the cylinder (green sectors). For iiq^' < vq < v^-^^ it neither 
reaches the bottom nor the piston (red sector), whereas for 

(2) (3) 

Vq ' < vo < Vq ' it reaches the piston but not the bottom 
(blue sector). 



previous velocities Vi,i = 0, .., {j — 1) which makes fur- 
ther analytical progress hard [12]. To disentangle this 
recursive dependence we introduce the logarithmic time 
variable 



T{t) = In 1 + t 



(4) 



For an expanding piston, the logarithmic time increases 
with real time, while it decreases for a compression. The 
logarithmic time increments between collisions of the par- 
ticle at the bottom are then given by 



At, ^ 



T{tj+i) ~T{tj) ^In 



2u 



for Vj < 2u, 
otherwise. 



(5) 

For Vj < 2u, the velocity of the particle after the colli- 
sion with the piston is 2u — Vj > such that the particle 
does not return to the bottom at all. We therefore set 
Atj = oo for this case. Since the velocities Vj after suc- 
cessive collisions with the bottom are independent identi- 
cally distributed random variables, this holds true for the 
time increments Atj as well. This fact will allow us to 
use methods from the theory of continuous time random 
walks [m [52] to determine the desired work distribution. 

In order to find the pdf p{W; r) that during a process 
of duration r = In L j / Lq the work W is performed we 
have to distinguish three types of particle trajectories (cf. 
Fig. [I] for the case of an expansion): 

1. The "null" case in which the particle neither hits 
the piston nor the bottom. This happens with 
probability Pq (t) that the initial velocity is between 
vj^'' — —xq/T and v^^ = {Lq — xq)/T + u. Since 
no work is performed in this case we have for the 
corresponding pdf pn(VF;r) = S{W)Po{t). 
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2. The "piston" case in which the particle reaches the 
piston and performs work but does not hit the bot- 

(2) 

torn. It occurs if the initial velocity lies between Vq ' 

and Uq — {2Lo — X(i)/T + 2u. The corresponding 
pdf will be denoted by Pp(W^; r). 

3. The "bottom" case in which the particle hits the 
bottom at least once and therefore gets thermal- 
ized. This happens if the initial velocity vq is either 
smaller than v^'^ or larger than v'"^\ We denote the 
pdf for this case by Ph{W] r). 

The total probability density that work W is per- 
formed in the time interval t is therefore given by 



p{W-T) = 5{W)Pa{T)+p^{W-T)+p^{W-T) . (6) 

The simplest case to analyze is the first one. From ^ 
we find 



r (2) 

Po{t) ^ / dxo / dwoPo(a;o,wo) 

Jo Ji;'^' 



1 



1 



A erf (A) + -^e-^^ - Serf (B) 



= e-" +Cerf(C) + — e 



(7) 



with the constants 



A 



-u, B 



2 Lo(e- - 1) 



(3 u 



2 ~ 1 



(8) 

In the second case the work W = —2u[vq — u) is per- 



formed. We hence get 




PpiW;T)^ dxQ dvo po{xq,vo)S[W + 2u{va~u)] 



W 
2u 
W 



exp 



'13 [u 



(e^ - 1) 1 



2v? 



2,1 



W 



(9) 



The most complex case is the last one. It consists itself 
of three ingredients: the initial period up to the first 
collision with the bottom, the middle period of repeated 
collisions with the bottom, and the final period following 
the last collision with the bottom. 

Let us start with the determination of the joint pdf 
Pi (A Wo, To) for the duration and work of the first period. 
The logarithmic time tq of the first contact of the particle 
with the bottom is given by 



To 



Lqvo 



In 1 



(2Lo-xo)u 
(vo-2u)Lo 



for va < min (o, ^) , 
for vq > max (2u, 



(10) 

The values of vq not covered above are those realiza- 
tions considered by Po(''") ^nd Pp{W; t). In the first case 
of ( 10 1 the particle moves directly to the bottom and 



hence does not transfer work before the first contact with 
the bottom. In the second case it reaches the piston first 
and does transfer work before hitting the bottom. From 
^ we get 



Pi{AWo,To) =5{AWo) 



dvop{xo,vo)5 



To 



In 1 



xqU \ 
Lovoj 



dxo 



e^o f S{AWo) 
'\u\V2^\ 



dvQ p{xo, vq)S[Wq + 2u{vo - u)] S 

2' 



To 



In 1 



Wo < nun 1^0 

^ (2Lo - xq)u 
(wo - 2it)Lo 



Xl^U 



Vq > max 2u, 



XqU 



2 



AWq 



2v? 



1 



exp 



exp 
2 



u 



e7« - 1 

2" 



AWo 
2u 



1[0<x(AWo,to)<1]^1[to>0] 



(11) 



r 



with x(AW^o,ro) = 2 + (e-o _ i)[MVa/{2u^) + i]. 

Note that / dAWo dToPi(AVFo, to) < 1 since the par- 
ticle never reaches the bottom if its initial velocity is 
between zero and 2u. This case is covered by Pq{t) and 
Pp{W\t) respectively. 



To characterize the middle period of repeated collisions 
we introduce the pdf ^{AW, At) giving the joint proba- 
bility density that between two collisions of the particle 
with the bottom the logarithmic time increases by At 
and the work AW is performed. It is most conveniently 
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specified by its Fourier transform in W and (bilateral) 
Laplace transform in At 



ij){k,X)= dAr / dAWe-^^^+'''^^il;{AW,AT). 

J — oo J — oo 

. (12) 

Since r is either positive for expansions and negative for 
compression we have to use the bilateral Laplace trans- 
form in contrast to the usual formalism which uses an 
unilateral Laplace transform [22l. The Laplace trans- 
form is defined for 5ftA > for u > and 5RA < for 
M < 0. Using ([T]), ([5]), and ([2| it acquires the explicit 
form 



max(2u,0) 



The integration starts firom max(2u, 0) since At{v) = 
00 for V < 2u (Eq. ([5|). For u > we have 

limA^oo V'lO, A) < 1 which is the probability that the 
particle returns to the bottom. 



The period after the last collision with the bottom is 
characterized by the fact that the time Arf remaining 
till the end of the process is not sufficient for the par- 
ticle to reach the bottom again. To determine the final 
contribution to the work we have to distinguish whether 
the particle reaches the piston a last time or not. If 
Arf < hi[v /{v — u)] it does not and there is no final work 
increment (first term in (14)). If \ti[v/{v — u)] < At{ < 
hi[v / {v — 2u)] the particle hits the pistons again and per- 



(13) 



forms a final contribution AWf to the work (second term 
in (14)). Note that Arf cannot exceed \n[v/{v — 2u)] ac- 
cording to ([5]). 



/"OO 

Pf (AWf ; Arf) =d{AW[) / dv (f>{v) 1 
Jo 



V < 



1 - e-^-'f 
+ I dv(f>{v)5[AWf + 2u{v-u)]t 



--5{AWi) <^ 1 - exp 



1 V 1 




<V < 



2u 



1 - e-^^f 



2uV\ 2u 



AWf 

<u — < 

2u 



2u 



(14) 



Now, assume that we start for r = at the bottom and 
denote the pdf for this process by Ps{W;t). When the 
particle hits the bottom again, the process can be de- 
scribed by the same distribution Ps{W; r) from this point 
on (also called a renewal) [21] . This leads to the equation 



p,iW;T)^MW;T) 



dW 



dr' V'(W, t')p^{W ~W'\t- r'). 

J 

(15) 



This relation is quite intuitive: the first part describes 
the case that we do not hit the bottom again, the second 
part describes the renewal property which is a process 
of remaining length r — r' which already has acquired a 
work contribution of W' in advance. Taking the Fourier 
transform with respect to W and the Laplace transform 



with respect to r turns Eq. (15) into the algebraic equa- 
tion 



which gives 



Ps(fc; A) 



Pf(fc; A) 
l-V'(fc,A)' 



(17) 



We are now in a position to derive an equation for the 
overall Pb(M^; t) characterizing the bottom case. We have 
to accommodate the fact that we do not necessarily start 
from the bottom, but the first hit is described by the pdf 
Pi(W, r). Along the lines of the renewal argument, we get 
the convolution 



Pb(W;r)= / dW \ At' p;{yV' ,T')p^{yV-W'\T-T'). 

J —00 J —00 

(18) 

Taking again the Fourier transform with respect to W 
and the Laplace transform with respect to r gives 



Pb(fc; A) 



Pi(fc, A)pf(fc; A) 



(19) 



Uk; A) = ^f (fc; A) + i,{k, X)Psik; A) (16) 



Upon inverse Fourier and Laplace transformation we get 
the result for p\^{W;t) which completes the analysis of 
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the bottom case and fixes the last term in ([6|. Unfortu- 
nately neither the u-integral in (131 nor the inverse trans- 
formation of (19 1 can be performed analytically. How- 



^'^^ with the free-energy change 



ever, they can be done efficiently using numerical meth- 
ods. 

In Figs. [2] and |3] we show several examples of work 
distributions obtained from ^ [inserting Eqs. ([T]), (|9]), 
([19]), ([TT]), ^ and ([14|] and compare them with results 
from numerical simulations of the particle dynamics. The 
agreement between analytical and numerical results is 
always very good. The relative deviations are typically 
below 1%. 

Fig. [2] shows the evolution of p{W; r) with time. For 
both compression and expansion the probability that the 
particle has hit the piston more than once is small for 
small times T. The distribution of work therefore shows 
only one major peak. With increasing T, also the proba- 
bility for multiple collisions increases, resulting in a dis- 
tribution with several peaks. At very long times, these 
peaks begin to broaden and finally merge. 

The distributions for expansion and compression 
mainly differ in the time evolution of the peaks at small 
work values. In an expansion the probability for the par- 
ticle to hit the piston only once does not decrease with 
time and hence the height of the peak in p{W, r) at small 
W remains roughly constant, cf. the curves for T = 20 
and T = 100 in Fig.[2](b). For a compression the steadily 
decreasing volume reduces the probability for just a single 
collision between particle and piston and correspondingly 
the peak at small values of W gradually disappears, as 
shown in Fig. [2] (a). 

The influence of the piston speed u on the work dis- 
tributions is shown in Fig. |3] For a very slow piston 
there are many collisions with the particle, the work in- 
crements of the individual collisions are independent ran- 
dom variables and by the central limit theorem p{W,t) 
becomes Gaussian. For intermediate piston speeds, the 
distribution develops multiple peaks originating from re- 
alizations with different numbers of collisions. For high 
piston speed, p{W, r) is dominated by a peak at = 
(not shown in the Figs.) stemming from realizations 
where the particle does not reach the piston at all and 
a broad unimodal contribution originating from just a 
single collision between particle and piston. 

The main difference between the distributions shown 
in Figs. |2] and [3] and those obtained for adiabatic com- 
pression and expansion respectively [9] lies in the long- 
time behavior. Whereas in the adiabatic case the fea- 
tures resulting from individual collisions leave their trace 
in p{W; r) for arbitrary long times they are washed out 
in the isothermal case as time goes by. This is consistent 
with the fact that the only random variables in the for- 
mer case are those from the initial condition whereas in 
the isothermal process a new random variable enters the 
stage after each collision with the bottom. 

To check our results for p{W; r) we have tested them 
against the Jarzynski identity [13] and the Crooks rela- 
tion [TS]. The Jarzynski equation stipulates {e~^^) = 



AF 



13 Lo P 



(20) 



Fig. [4] demonstrates that the average {e~^^) determined 
with various distributions p{W; r) obtained from ([6| in- 
deed coincides with e~^^^ of the corresponding process. 
The deviations are less than 1% and can be assigned to 
round-off errors in the numerics for the negative tails of 
piW;T). 

The Crooks relation 



p{-W;t) 



(21) 



involves the distribution of work for the reverse process, 
p{W;t). The reverse process of an expansion (compres- 
sion) from Lq to L{ at piston speed u is a compression 
(expansion) form Lf to Lq at piston speed —u. In Fig. [s] 
the logarithmic ratio l3~^\n{p{W;T)/p{—W;T)) is plot- 
ted against W. Deviations from ( [2l| ) are only visible 
where p{W; t) or p{—W\ t) are very small (less than 
10^^) and can again be attributed to round-off errors. 



B. Quasi-static limit 

While it seems not possible to get an explicit analytic 
expression for Eq. ([6]), it is feasible to get results in the 
limit u — ^ 0. In order to shorten the derivation we make 
the further assumption that the particle starts at the 
bottom of the cylinder. The error resulting for the work 
is of order u and can be neglected. 

For this process the probability Ps{W]t) can be ex- 
pressed in Fourier-Laplace space as (Eq. (17)) 



Ps(fc; A) 

With the definition 
9u{k,X) 




Pf{k;X) 
l-V'(fc,A)' 

1 j cxY){—ikuv) 



(22) 



(23) 



A — iu. 



-k + o{u) 



we have from Eq. (|13|) 
i'{k,X) 



2iku' 



g2u{k, A) 



(24) 



and from Eq. ( 14 ) 



Pf{k, A) = ^ (1 - .g„(0, A) + e2^'=^'ff„(2fc. A) - ^(fc. A) 



uy/2'Kl3 + o{u). 



(25) 
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Figure 2: Time evolution of the distribution of work p{W; r) for different durations T — Lo(e^ — 1)/^) for botfi compressions (a) 
and expansions (b). Tlie system starts in equilibrium at temperature l/P — 1, the initial length of the cylinder is set Lq = 1, 
the piston speed is ti = ±0.2. Shown are the analytic expressions [lines, eq. ^] together with results of molecular dynamics 
simulations (circles) of 10^ trajectories each. The (5-peaks at = are not displayed. 




Figure 3: Distribution of work p{W; r) for compressions (a) and expansions (b) starting from equilibrium at temperature 
= 1 for different speeds u of the piston. The length of the cylinder is changed from Lq = 1 to Lf — 2 and vice versa 
respectively. Shown are the analytic expressions [lines, eq. (j6|] together with results of molecular dynamics simulations (circles) 
of 10^ trajectories each. The 5-peaks at Vl^ = are not displayed. 



Plugging this into Eq. ( 22 ) results in 



collisions stay correlated and the work fluctuates even in 
the quasi-static limit. 



(26) 



The inverse Fourier-Laplace-transform gives the final 
result 



Ps{W;T) = S{W + r'r) 
^6{W -AF). 



(27) 



This result is calculated for the simplified process 
Ps{W; r), however the corrections to p{W; r) is only the 
behavior till the first hit at the bottom which is of lower 
order than considered here. Hence, in the quasi-static 
limit and to leading order, the work does not fluctu- 
ate and is always equal to the free-energy change AF. 
This result is also suggested by the law of large numbers 
since the work is performed in infinitely many collision 
between the particle and the piston. In contrast, during 
an adiabatic compression or expansion [3] the individual 



IV. THE FULL CYCLE 



We now combine our results of section ITlI Al for isother- 
mal volume changes with those of Lua and Grosberg [S] 
for adiabatic ones to analyze the fluctuations of work and 
heat in a full cyclic process consisting of four strokes, sim- 
ilar to a Carnot process. The organization of the cycle 
together with the relevant notations is shown in Fig. [6] 
We choose the cycle to start with an isothermal com- 
pression at inverse temperature /3coid- It is followed by 
an adiabatic compression changing the inverse tempera- 
ture from /3coid to /3hot- The third stroke is an isothermal 
expansion at the higher temperature followed by an adi- 
abatic expansion back to initial volume. This definition 
of the cycle is somewhat arbitrary; the cycle could start 
at any stroke. As before the initial position and velocity 
of the particle are denoted by xq and vq respectively, Xi 
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Figure 4: Verification of the Jarzynski equality. The mean 
{e~^^) for compressions and expansions starting from Lo = 1 
in equifibrium at temperature 1 is determined for different 
speeds of the piston and different final lengths of the cylinder. 
The results are plotted against e~^^^ which is proportional 
to the final length L{. 




Figure 6: Variables used to describe the full cyclic process 
discussed in section |IV[ Xi and Vi denote the position and 
the speed of the particle at the end of the i-th stroke, Li is 
the length of the cylinder at this moment, and Wi and Qi 
denote the work and heat respectively exchanged during the 
i-th stroke. 




of the different strokes are not known. Note that it is not 
even acceptable to use ^ for the initial distribution of 
(xq, vq)^ since in general we would not come back to it af- 
ter the fourth stroke, and would therefore not describe a 
time-periodic steady state (TPSS). We hence start with 
an investigation of the quasi-static limit of the cycle in 
which the situation is somewhat simpler and some com- 
parison between analytical work and simulations is still 
possible. After this we elucidate the general case with 
arbitrary values of Uc and on the basis of numerical 
simulations. 



Figure 5: Verification of the Crooks relation. The log-ratio 
of the work distribution p{W) of an isothermal compression 
from Lo = 2 to Lf = 1 and p{—W) of the corresponding 
reverse process is plotted for different values of u. All points 
lie on the line W — IS.F as required by (21 1. For |m| = 0.02 and 
W > 3, p(W) and p(—W) are less than 10~^® and dominated 
by round-off errors. Hence these values are not shown. 



and Vi are their respective values after the i-th stroke. 
The length of the cylinder at these times is denoted by 
Li. To keep the analysis simple we choose the same pis- 
ton speed Uc for both compression strokes and the same 
value Uc for both expansions. 

The heat transferred during the isothermal strokes can 
be determined by applying the first law of thermodynam- 
ics to each individual realization of the process. Since the 
change in internal energy is just the change in kinetic en- 
ergy of the particle, we find 



^cold 



Q 



hot 



1 ( 2 

2 ^ 1 

1 I 2 

2 (-3 



(28) 
(29) 



A general analysis of the full cycle is hampered by the 
fact that the distributions of Xi and Vi at the beginnings 



A. Quasi-static case 

The quasi-static case is defined by the combined limit 
Mc — > and Me — > 0. From section [ill B| we know that in 
these limits the work performed in the isothermal strokes 
equals the free-energy difference, cf. (27), i.e. 



lid 



1 ^^^1 

/3cold Lq 



= AFhot = --^ln^ 

Phot J^2 



(30) 
(31) 



In the adiabatic strokes the number of collisions Ui is 
given for small piston speed u by [9] 



2u y L,+i 



0(1). 



(32) 



After these collisions the particle speed is reduced to 

, Li 



2n,u 



' L 



+ 0{u) (33) 
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Figure 7: Comparison of the analytical solutions for the quasi-static case with simulations of 10^ cycles with a piston speed 
u = ±10^^. In the subplots, the respective joint pdfs /icoid(H^, Qcoid) (a) and /ihot(VK, Qhot) (b) of the work W and the heat 
Qi are shown in the center. The analytical solution is depicted by lines of equal probability, while each simulated cycle is 
represented by one black dot. The marginal probability pw{W) is shown on the top of each subplot whereas the marginal 
probabilities qi{Qi) can be seen at the right hand sides of (a) and (b). Here lines denote the analytical solutions, while the 
simulations are depicted as circles. The inverse temperatures are P^^^^ = 2 and I3^J^^ = 1, and in the isothermal stroke at the 
higher temperature the cylinder length changes from L2 = 1 to L3 = 2. 



and the work performed is therefore given by 
Wi+i ^ -2ntu{\vi\ - Uiu) 



2 \L^+, 



1 +0{u) 



(34) 



From ( 33 ) it follows that if Vi obeys a Maxwell distribu- 
tion at inverse temperature Pi then w^+i will be Maxwell 
distributed as well, however at inverse temperature 



•i+l 



(35) 



as expected for a classical ideal gas in one dimension. 
Since Q ensures that the Maxwell distribution is repro- 
duced after each isothermal collision we find back the 
well-known fact that in the quasi-static limit the system 
is always in equilibrium. 



Hence 




and we have to leading order 



W2 = ^ 



Wa 



/Shot 
/3hot 



1 



'cold 



(36) 



(37) 



(38) 



From ([281, Q, ([31]), p7|), and Q it is clear 

that the total work W — ^ '^^ll Qcoid 

and Qhot depend on the three velocities vq, wi, and 
only. These velocities are independent samples from the 
Maxwell distribution. The velocity V2 is then fixed by 

«l\//3cold/^hot- 



V2 

Thepdfp(W^,Qcoid , Qhot ) 
cle is given by 



characterizing the whole cy- 
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, ^fcold, Whot) 



/5coldV/3hot 

(27r)3/2 
<5 



dvi e 



dv3 e 



-/3ho 



'cold 



/3hot 

X 5 { Qcoid 



oo 
„,2 



-II- Ai^^hot 

,2 



Jhot 



- 1 



vt - V, 







2 V /^cold 

A-Fcoid (5 Qhot 



,,2 ffcold 

-'l /3hot 



AF, 



hot 



(39) 



To disentangle this expression we introduce the char- 
acteristic function 

p{k, k', k") = (exp {ikW + ik'Q,,id + «fc"Qhot}> (40) 

to get 



p{k,k',k") 



1 



^l-P'-hk~P~l^k' + P~^\^k" 



exp In (^) - /3coM^' + /3hotfc") } 



with = /3: 



hot- 



(41) 



The inverse Fourier- 



transformation of (411 cannot be done exactly in full 
generality. However, we may get analytical results for 
some marginal distributions. 

To begin with we consider the distributions of heat 
transferred in the individual strokes which we denote by 
<?coid((5coid) and ghot((9hot) respectively. Their charac- 
teristic functions are given by gcoid(fc) — p(0,fc,0) and 
9hot(fc) = p(OiO, fc) respectively. The inverse Fourier- 
transform of these characteristic functions can be done 
analytically and gives 

q,{Q^) = ^Ko [A(Q. + A^^,)] , i = cold, hot, (42) 

where Kg denotes a modified Bessel function of the sec- 
ond kind. A similar expression was found for the equi- 
librium heat fluctuations of a Brownian particle in a 
parabolic potential [531 123 . Related is also the result of 
an exponential tail of the distribution of injected power 
for a Brownian particle in contact with two heat reser- 
voirs as derived in 28 . In [27j it was shown that the 
distribution of heat has exponential tails even for an over- 
damped Brownian particle in a moving parabolic poten- 
tial, i.e. in a non-equilibrium situation. 

For our system also the joint pdf ^hot{W, Qhot) for the 
total work and heat exchanged with the cold reservoir 
can be determined exactly. Its characteristic function 



is given by /ihot(fc,fc') — p(fc, 0,fc') which upon inverse 
Fourier transformation yields 

Mhot(W^, Qhot) = '?hot((3hot) S{W + rjcQhot)- (43) 

Here ijc = I — Phot / /^coid is the Carnot efficiency. Inte- 
grating over Qhot we also find the marginal distribution 
Pw{W). Finally, /icoid(M^, Qcoid) is obtained by perform- 
ing the inverse Fourier transform numerically. 

In Fig. [7]/Zcoid(T^, Qcoid) and Athot(W, <9hot) are com- 
pared with results from molecular dynamics simulations. 
There is good agreement, with small deviations due to 
the fact, that the simulations were performed with fi- 
nite piston speed u = 10~^. As can be seen, also in the 
quasi-static limit, in a Carnot cycle with a single parti- 
cle, there remain substantial fluctuations of W, Qhot, and 
Qcoid- This seems surprising since both work and heat 
are transferred in infinitely many collisions of the parti- 
cle with the piston and the bottom, respectively. In fact, 
in the isothermal strokes the work is equal to the differ- 
ence in free energy and therefore does not fluctuate. In 
the adiabatic strokes, however, the individual collisions 
with the piston are all correlated with each other, the 
law of large numbers does not apply, and fluctuations in 
the work remain. It is these fluctuations from the adia- 
batic strokes that entail the non-trivial structure of the 
distribution for the total work. 

The heats exchanged with the reservoirs during the 
isothermal strokes are equal to the difference between the 
kinetic energy of the particle and the work performed, 
cf. (28) and ((29|). In the quasi-static limit the work in 



the isothermal strokes does not fluctuate and the heats 
therefore just depend on the two random velocities of the 
particle at the beginning and at the end of the respec- 
tive stroke. The influence of these velocities on the heats 
is not diminished by m — >■ and consequently the cen- 
tral limit theorem is again not applicable. Therefore, the 
distributions of Qhot and Qcoid are non-trivial as well. 
If there are many particles in the cylinder their inde- 
pendent contributions combine and the results of macro- 
scopic thermodynamics are, of course, restored. 

Also the correlations between heat and work are inter- 
esting. Since W and Qhot depend on the same random 
variables, namely vi and V2, they are proportional to each 
other as expressed by (43). Note that this relation also 



implies that the fluctuation theorem for heat engines as 
derived in [23 is trivially fulfilled.. 
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In contrast, Qcoid depends in addition on vq and the 
correlation with W is not as strong as for Qhot- We need 
to emphasize, however, that there is some arbitrariness 
to these statements since they depend on where we de- 
fine the beginning of our cycle. So far we have chosen the 
isothermal compression at the lower temperature as the 
first stroke, cf. Fig. [6] If, instead, we were starting with 
the isothermal expansion W and Qcoid were proportional 
to each other and W and Qhot would show only weak cor- 
relations. This is also the general picture for the quasi- 
static case: irrespective of the starting point of the cycle 
there is always one heat that is proportional to the work 
and one that shows relative fluctuations. This depen- 
dence on the precise definition of the cycle is restricted 
to the correlations, however, the marginal distributions 
of work and heat remain unchanged. 



B. Arbitrary piston speed and efficiency at 
maximum power 

We now come to the general case of arbitrary Me and 
Uc, where the TPSS operation is of particular inter- 
est. By definition, for the TPSS operation the relation 
p{xi,vi) = p{x4, V4) must hold. It seems to be difficult to 
determine these distributions analytically. We therefore 
rely on numerical results to characterize the TPSS. To 
this end we start simulations in thermal equilibrium and 
let them run, until the distributions pvj{W), 9hot(<5hot) 
and qcoid(Qcoid) no longer change. All numerical data 
discussed in this subsection were obtained after such a 
transient phase. Results for the distributions of work 
and heat are shown in Fig. [8] for three different speeds of 
the piston. 

In contrast to the quasi-static limit, for finite and 
Uc, the relations Wi = AFcoid and W3 = AF^ot no longer 
hold. Yet, if the piston speed is not too large the work is 
still transferred in many independent collisions the num- 
ber of which decreases with increasing Uc and Uc- By 
virtue of the central limit theorem one would hence ex- 
pect that for small piston speeds the work variables Wi 
and 11^3 become Gaussian distributed with their width 
increasing with increasing piston speed. This broadening 
of the distributions pw{W), QcoidlQcoid), and ghot((3hot) 
can be clearly seen by comparing Fig. [7] with Figs.[8f^a-d). 
When the piston speed increases further the probability 
for the particle to neither reach piston nor bottom dur- 
ing an entire cycle becomes noticeable. These realizations 
do not transfer work or heat and give rise to J-peaks in 
the distributions of W, Qcoid and Qhot, cf. Fig. [sjjc-f). 
Moreover, typical particle trajectories now involve a few 
collisions only and the distributions of work and heat ex- 
hibit a much more spiky structure than those from the 
case of slow driving, see Figs.[8](e) and (f). Nevertheless 
we again find exponential tails for the marginal heat dis- 
tributions shown in Figs, [s] similar to the results of [27] 
for a Brownian particle in a moving parabolic potential. 
We can no longer test our results against the fiuctuation 



theorem for heat engines f5S^ since the latter is valid only 
for cycles starting in equilibrium. 

Based on the average power (W) delivered by the cycle 
and the average heat fiuxes (Qcoid) and (Qhot) we may 
distinguish three different regimes of operation as shown 
in Fig. [9] The mean power (W) is defined as the total 
average work per cycle (W) divided by the duration of 
one cycle and the heat fluxes (Qcoid) and (Qhot) as the 
average total heat transfers per cycle divided by the du- 
ration. In (A) the piston moves slowly and the engine 
produces power. The efficiency 77 ~ — (W^)/(Qhot) starts 
from the Carnot value 77c in the quasi-static limit and 
decreases monotonically down to zero at the border of 
this regime where (W) =0. In the intermediate range 
(B) of piston speeds, the engine consumes work and heat 
is transferred from the hot to the cold bath. Here rj is 
negative. If the piston moves still faster, (C), work is con- 
sumed and both reservoirs are heated up. In this regime, 
7] is formally larger than one, but the device is clearly 
useless as a heat engine. 

Fig. [9] also shows, that the power delivered by the en- 
gine has a maximum at an intermediate value of Uc — 
—Uc- We have determined the efficiency at maximum 
power (EMP) rj* of our model engine for nine different 
combinations of /3hot and /3coid- For each combination we 
varied u^. and —Uc independently between 0.01 and 0.35 
to find the point of maximum power. The results are 
shown in Fig. [lO] together with the Curzon-Ahlborn effi- 
ciency r/cA = 1 — \/l ~ Vc [13] and the bounds rj- — r]c/2 
and 77-1. = T]c / {"2 — r]c) derived in Schmiedl [19j, Espos- 
ito ct al. [20]. Our findings are consistent with these 
bounds and almost saturate the upper one. The devia- 
tion from the Curzon-Ahlborn efficiency is not surprising 
since these authors assumed a well defined temperature 
of the working medium whereas our simulations show 
that the velocity distribution p{v) at maximum power is 
markedly different from a Maxwell distribution. We also 
note that the numerical results for the EMP reported in 
[T6] are somewhat larger than ours thereby violating the 
upper bound rj^. This might be due to the fact that in 
this study the optimization of the power was restricted 
to Mo = — Mc, a relation that our optimized values for Ue 
and Uc do not fulfill. 



V. SUMMARY 

In the present paper we gave a detailed analysis of 
the stochastic energetics of isothermal compressions and 
expansions of a classical ideal gas consisting of a single 
particle. In our model the volume of a one-dimensional 
cylinder is changed by a piston moving with constant 
speed u. The enclosed particle performs work by elastic 
collisions with this piston and exchanges heat in inelastic 
collisions with the bottom of the cylinder. The kinetic 
energy of the particle is of the order of the average ther- 
mal energy per degree of freedom. Accordingly, both 
work and heat are strongly fluctuating variables, and our 
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Figure 8: Distribution of work W and heat Qcow and Qhot exchanged with the cold and the hot bath respectively during 
one cycle for different piston speeds Uo ~ — Uc- In each subplot the center plot displays the joint probability ^coid(Wi Qcoid) 
or fJ-hotiW, Qhot) respectively. The marginal work distribution pw{W) is always shown on the top whereas the marginal 
distributions gcoid(Qcoid) and (;hot(Qhot) are given at the right. Every black dot represents one of 10^ simulated cycles. Parameter 
values are /3^Jj = 1, /^^ot ~ ^nd Me = — Mc = 0.01 (a,b), Uc — —Uc = 0.1 (c,d), and Uc — — Mc ~ 1 (e,f). In all cases the length 
of the cylinder changes in the isothermal stroke at the higher temperature from L2 = 1 to L3 = 2. The solid line in subfigures 
b,d, and f indicates the quasi-static limit. 



focus was on the determination of their probability distri- 
butions. The piston speed is arbitrary, and in addition to 
quasi-static volume changes we were in particular inter- 
ested in situations where u is comparable to the typical 
velocity of the particle. 

We first analyzed isothermal compressions and expan- 
sions. Introducing a logarithmic time scale we were able 
to disentangle correlations between successive collisions 
of the particle which had hampered analytical progress 
so far. Using elements of renewal theory we then calcu- 
lated the characteristic function of the work distribution 
analytically. The distribution itself was determined using 
a numerical implementation of the inverse Fourier trans- 
form. The results agree very well with molecular dynam- 
ics simulations and fulfill the Jarzynski equality as well as 
the Crooks relation. In the quasi-static limit we recover 
the result that the work ceases to fluctuate and coincides 



with the difference in free energy between the final and 
the initial state. For increasing piston speed we first find 
a broadening of this (5-distribution into a Gaussian before 
at even larger values of u the typical number of collisions 
becomes of order one and a multi-modal shape of the 
work distribution emerges. 

We then combined these findings with the results of 
Lua and Grosberg for adiabatic volume changes to ana- 
lyze a Garnot-like cycle involving two isothermal strokes 
at different temperatures and two adiabatic strokes con- 
necting them. For the quasi-static limit we derived an- 
alytic results for the joint distribution of the total work 
and the heat exchanged with the hot reservoir, as well 
as for the characteristic function of the heat exchanged 
with the cold bath. Again the results are in very good 
agreement with molecular dynamics simulations. Some- 
what surprisingly strong fiuctuations in work and heat 
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Figure 9: (a): Mean total power (W) and mean heat fluxes (Qcoid.hot) during the hot and the cold isothermal stroke respectively 
in dependence of the piston speed u. The simulations are conducted for 2- 10^ cycles each, with I3~^^ = 1 and fi^^^ — 2, and the 
position of the piston changing from L2 = 1 to L3 = 2 during the isothermal expansion at /3j^Jj. (b): Average efficiency obtained 
from the same simulations. There are three regimes of operation: {W} < (A), (W) > and (Qhot) > (B), (Qhot) < (C). 
(c): Regimes of operations as obtained from simulations of 10^ combinations of piston speeds Uc and Uc- 




Carnot efficiency rjc 



Figure 10: Comparison of the efficiency at maximum power 
rj* (red circles) for different temperature combinations with 
the Curzon-Ahlborn-efficiency rjcA (solid line) and the lower 
and upper bounds rj- (dashed-dotted line) and 77+ (dashed 
line) obtained in [191 I20j . The power was optimized using 
simulations of 10^ cycles each for which the cylinder length 
again changes from L2 = 1 to I/3 = 2 during the isothermal 
expansion at the higher temperature. 

remain even in the quasi-static limit: only when consid- 
ering a large ensemble of independent particles the results 
of macroscopic thermodynamics are recovered. 

Of special interest are again situations beyond the 
quasi-static regime. Unfortunately, it seems very hard to 
determine the invariant density for particle position and 



velocity of the time periodic steady state describing the 
stationary operation of the cycle. Therefore we had to 
rely exclusively on numerical simulations to explore the 
non-equilibrium performance of the cycle. For increas- 
ing piston speed the strong correlations between work 
and heat observed in the quasi-static limit get weaker. 
The work delivered per time increases and at intermedi- 
ate values of the piston speed there is a maximum of the 
power output. We determined the efficiency at maximum 
power for several values of the parameters and showed 
that they are consistent with recently derived bounds. 
When the piston speed increases further a complicated, 
spiky shape of the probability distributions emerges. At 
the same time the efficiency deteriorates until finally the 
cycle ceases to work as a heat engine altogether: on aver- 
age it then consumes work and heats up both reservoirs. 

It should be interesting to extend our investigations 
to the case where the position of the piston is not de- 
terministically prescribed but is itself a fluctuating vari- 
able. It would then be possible to analyze the influence 
of intermingled forward and backwards steps which are 
characteristic for nanoscopic machines. 



Acknowledgments 

We would like to thank Chris Van den Broeck and 
David Lacoste for stimulating discussions. 



[1] U. Seifert, Rep. Prog. Phys. 75, 126001 (2012). 



[2] C. Jarzynski, Annu. Rev. Condens. Matter Phys. 2, 329 



13 



(2011). [17 
[3] M. Esposito, Phys. Rev. E 85, 041125 (2012). 
[4] U. Soifcrt, Phys. Rev. Lett. 106, 020601 (2011). [18 
[5] C. Van den Broeck, N. Kumar, and K. Lindenberg, 

Phys. Rev. Lett. 108, 210602 (2012). [19 
[6] V. Blickle and C. Bechinger, Nat. Phys. 8, 143 (2011). 
[7] O. Mazonka and C. Jarzynski, [20 
arXiv:cond-mat/9912121vl [cond-mat.stat-mech] 
(1999). [21 
[8] R. van Zon and E. G. D. Cohen, Phys. Rev. E 67, 

046102 (2003). [22 
[9] R. C. Lua and A. Y. Grosberg, J. Phys. Chem. B 109, 

6805 (2005). [23 
[10] L Bena, C. Van den Broeck, and R. Kawai, Europhys. 

Lett. 71, 879 (2005). [24 
[11] B. Cleuren, C. Van den Broeck, and R. Kawai, Phys. 

Rev. E 74, 021117 (2006). [25 
[12] A. Baule, R.M.L. Evans, and P.D. Olmsted, Phys. Rev. 

E 74, 061117 (2006). [26 
[13] R. Noltc and A. Engcl, Physica A 388, 3752 (2009). 
[14] C. Jarzynski, Phys. Rev. Lett. 78, 2690 (1997). [27 
[15] G. E. Crooks, Phys. Rev. E 61, 2361 (2000). 
[16] Y. Izumida and K. Okuda, Europhys. Lett. 83, 60003 [28 
(2008). 



Y. Izumida and K. Okuda, Phys. Rev. E 80, 021121 

(2009). 

M. Esposito, K. Lindenberg, and C. Van den Broeck, 

Phys. Rev. Lett. 102, 130602 (2009). 

U. Schmiedl, T. Seifert, Europhys. Lett. 81, 20003 

(2008). 

M. Esposito, R. Kawai, K. Lindenberg, and C. Van den 
Broeck, Phys. Rev. Lett. 105, 150603 (2010). 

E. W. Montroll and G. H. Weiss, J. Math. Phys. 6, 167 
(1965). 

R. Balescu, Statistical Dynamics - Matter out of 
Equilibrium (Imperial College Press, London, 1997). 

F. L. Curzon and B. Ahlborn, Am. J. Phys. 43, 22 

(1975). 

D. Chatterjee and B. J. Cherayil, Phys. Rev. E 82, 

051104 (2010). 

A. Imparato, L. Peliti, G. Pesce, G. Rusciano, and 
A. Sasso, Phys. Rev. E 76, 050101 (2007). 
N. A. Sinitsyn (2011), arXiv:1111.7014vl 
[cond-mat.stat-mech] . 

R. van Zon and E. G. D. Cohen, Phys. Rev. E 69, 
056121 (2004). 

P. Visco, J. Stat. Mech. 2006, P06006 (2006). 



